Temporal Dynamics of Photon Pairs Generated by an Atomic Ensemble 
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The time dependence of nonclassical correlations is investigated for two fields (1,2) generated by 
an ensemble of cold Cesium atoms via the protocol of Duan et al. [Nature 414, 413 (2001)]. The 
correlation function R(t\,t2) for the ratio of cross to auto-correlations for the (1,2) fields at times 
(ti, £2) is found to have a maximum value JJ max = 292 ±57, which significantly violates the Cauchy- 
Schwarz inequality R < 1 for classical fields. Decoherence of quantum correlations is observed over 
Td — 175 ns, and is described by our model, as is a new scheme to mitigate this effect. 
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In recent years quantum measurement combined with 
conditional quantum evolution has emerged as a powerful 
paradigm for accomplishing diverse tasks in quantum in- 
formation science [j], Qj y, 01 . ^ or example, Duan, Lukin, 
Cirac and Zoller (DLCZ) |4| have proposed a scheme for 
the realization of scalable quantum communication net- 
works that relies upon entanglement created probabilisti- 
cally between remotely located atomic ensembles. By uti- 
lizing successful measurements to condition subsequent 
steps in their protocol, DLCZ have developed a scheme 
that has built-in quantum memory, entanglement purifi- 
cation and resilience to realistic sources of noise, thereby 
enabling a quantum repeater architecture to overcome 
photon attenuation 0, Q ■ 

Central to the DLCZ protocol is the ability to write 
and read collective spin excitations into and out of an 
atomic ensemble, with efficient conversion of discrete spin 
excitations to single-photon wavepackets. Observations 
of the resulting non-classical correlations between the op- 
tical fields generated from writing and reading such spin 
excitations have recently been reported by several groups, 
both at the single-photon level yj la, Lj as appropriate to 
the protocol of DLCZ and in a regime of large photon 
number n ~ 10 3 — 10 7 10]. Generation and detection 
efficiencies have now been improved so that excitation 
stored within an atomic ensemble can be employed as a 
controllable source for single photons Q. 

A critical aspect of such single-photon wavepackets is 
that they are emitted into well-defined spatio-temporal 
modes to enable quantum interference between emissions 
from separate ensembles (e.g., for entanglement based 
quantum cryptography 4i]). However, with the excep- 
tion of the verification of the time-delay implicit for 
the Raman processes employed [l(| , experiments to date 
have investigated neither the time or spatial dependence 
of quantum correlations for the emitted fields from the 
atomic ensemble. The high efficiencies achieved in Ref. 
J8J now enable such an investigation for the temporal 
properties of nonclassical correlations between emitted 
photon pairs, which we report in this Letter. 

Specifically, we study the time dependence of quantum 
correlations for photons emitted from an ensemble of cold 
Cesium atoms, with photon pairs created sequentially by 
classically controlled write and read pulses. Large vio- 
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FIG. 1: (a) Schematic of experiment. Write and read pulses 
propagate into a cloud of cold Cs atoms (MOT) at times Tj 
and Tj ± At respectively, and generate the correlated output 
fields 1 and 2. Quantum correlations for these fields at times 
(£1, £2) are investigated by way of photoelectric detection, (b) 
The relevant atomic level scheme. 



lations of the Cauchy-Schwartz inequality R < 1 for the 
ratio of cross to auto-correlations are observed for the 
generated fields, with i? max = 292 ±57 ^ 1. By contrast, 
previous measurements have reported violations R ^ 1 
only for detection events integrated over the entire dura- 
tions of the write and read pulses (R — 1.84±0.06 in Ref. 
0, R = 1.34 ± 0.05 in Ref. Q, and R = 53 ± 2 in Ref. 
|3). We also map the decay of quantum correlations by 
varying the time delay between the write and read pulses, 
and find a decoherence time Td — 175 ns. We have devel- 
oped a model to describe the decoherence and find good 
correspondence with our measurements. This model is 
utilized to analyze a new proposal that should extend 
the correlation times to beyond 10 fj,s, which would allow 
for entanglement between atomic ensembles on the scale 
of several kilometers. 

Our experimental procedure is illustrated in Figure ^ 
An initial write pulse at 852 nm creates a state of collec- 
tive excitation in an ensemble of cold Cs atoms, as her- 



aided by a photoelectric detection event from the Raman 
field 1. After a user-programmable delay At, a read pulse 
at 894 nm converts this atomic excitation into a field exci- 
tation of the Raman field 2. The (write, read) pulses are 
approximate coherent states with mean photon numbers 
(10 3 ,10 5 ), respectively, and are focussed into the MOT 
as TEMoo Gaussian beams with orthogonal polarizations 
and beam waist wq = 30 /im. This scheme is imple- 
mented in an optically thick sample of four-level atoms, 
cooled and trapped in a magneto-optical trap (MOT) 
|l l| . In particular, we utilize the ground hyperfine levels 
65i /2 ,F = {4; 3} of atomic Cs (labelled {[a); |6)}), and 
excited levels {6P 3/2 , F = 4; 6P 1/2 , F = 4} of the D 2 ,D 1 
lines at {852; 894} nm (labelled {|e); |e'}}). 

Each attempt to generate a correlated pair of photons 
in the (1, 2) fields is preceded by shutting off the trapping 
light for 700 ns. The re-pumping light is left on for an 
additional 300 ns in order to empty the \b) state, thus 
preparing the atoms in |o). The 7 th trial of a protocol is 
initiated at time Tj when a rectangular pulse from the 
write laser beam, 150 ns in duration (FWHM) and tuned 
10 MHz below the [a) — > |e) transition, induces sponta- 
neous Raman scattering to level |6) via \a) — ► |e) — > \b). 
The write pulse is sufficiently weak so that the probabil- 
ity to scatter one Raman photon into a forward propa- 
gating wavepacket is much less than unity for each pulse. 
Detection of one photon from field 1 results in a "spin" 
excitation to level \b), with this excitation distributed in 
a symmetrized, coherent manner throughout the sample 
of N atoms illuminated by the write beam |J . Regard- 
less of successful detection of a photon in field 1, we next 
address the atomic ensemble with a read pulse at a time 
Tj + At, where At is controlled by the user. The read 
light is a rectangular pulse, 120 ns in duration, tuned to 
resonance with the \b) — ► |e') transition. 

To investigate the photon statistics, we use four 
avalanche photodetectors, a pair for each field i, labelled 
as DiA,iB, which are activated at (Tj,Tj + At) with 
i = (1,2), respectively, for 200 ns for all experiments. 
The quantity p T (ti,t m ) is defined as the joint probabil- 
ity for photoelectric detection from field I in the interval 
[Tj + t; , Tj +ti+r] and for an event from field m in the 



interval [Tj + 1,, 



) -* j T £7; 



-t], where I and m equal 1 or 2. 



Pr{ti, t m ) is determined from the record of time-stamped 
detection events at D\a, is, F>2A.ib ■ In a similar fash- 
ion, q T (ti, t m ) gives the joint probability for detection for 
fields (I, m) in the intervals ([Tj + t u Tj + ti + r], [T k + 
t m 1 Tk + t m + t] ) for two trials k 7^ j. 

Following Refs. pj Q, we introduce the time- 
dependent ratio 7? T (ii,i 2 ) of cross-correlation to auto- 
correlation for the (1,2) fields, where 



«,<*,,*,). Mtl ' !j,li 



Pr(tl,ti)p T (t 2 ,t 2 ) 



(1) 



This ratio is constrained by the inequality R T (t\, t 2 ) ^J 1 
for all fields for which the Glauber-Sudarshan phase- 
space function is well-behaved (i.e., classical fields) 
j . Beyond enabling a characterization of the quan- 




FIG. 2: Probability for joint detection from the fields (1, 2) at 
times (ti,ta) with origin at the beginning of the write pulse, 
(a) p T (ti,t2), and (b) q T (t\,t2) for overlapped write and read 
pulses, At = 50 ns, with the solid line corresponding to t-z = 
£1. (c) p T (t\,t2) and (d) gr(ti,ta) for consecutive write and 
read pulses, At = 200 ns. In all cases, the bin size r = 4 ns, 
and the joint probabilities p T and q T have been scaled by 10 9 . 



turn character of the (1, 2) fields in a model independent 
fashion, measurements of i? T (ti,t 2 ) also allow inferences 
of the quantum state for collective excitations of single 
spins within the atomic ensemble. 

The first step in the determination of R T (t\, t 2 ) is the 
measurement of the joint probability p T (ti,t 2 ) for the 
(1,2) fields, and for comparison, q T (ti,ta) for indepen- 
dent trials. In our experiment, we focus on two cases: 
(/) nearly simultaneous application of write and read 
pulses with offset At = 50 ns less than the duration of 
either pulse, and (77) consecutive application of write 
and read pulses with At = 200 ns longer than the write, 
read durations. Results for p r (ti,t 2 ) and q T (ti,t 2 ) are 
presented in Fig. [3 as functions of the detection times 
(ti,t 2 ) for the fields (1,2). For both At = 50 and 200 
ns, p T (t\^2) ^> <7r(ti,t 2 ), indicating the strong corre- 
lation between fields 1 and 2, with the maximal ratio 
ff[, 2 (ti,t 2 ) = p T (ti,h)/q T (tiM) ^^2>|-Sf nicn is mucn 
greater than reported previously 0, H @- I n Fig- HI 
t = 4 ns, leading to statistical errors of about 8% for the 
largest values shown. 

In case (7) for nearly simultaneous irradiation with 
write and read pulses, Fig. [3(a) shows that p T (fi,t 2 ) 
peaks along the line t% — 1\ = 5ti2 — 50 ns with a width 
At 12 — 60 ns, in correspondence to the delay 8t\2 and 
duration Ati 2 for read-out associated with the transition 
16) — > \e') — ► \a) given an initial transition \a) — > |e) — » \b) 
jlOj . Apparently, the qualitative features of p T (ti,t 2 ) de- 
pend only upon the time difference between photon de- 
tections in fields 1 and 2 (i.e., p r (ti,t 2 ) « F(t 2 — £1)). 
In case (77) with the read pulse launched 200 ns after 
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FIG. 3: Probability of joint detection p T (ti,ti) for field I 
(squares) and p T {ti,t-i) for field 2 (diamonds) as functions of 
respective detection times £i and £2. Bin size r = 30 ns. (a) 
Overlapping write and read pulses, At = 50 ns; (b) Consecu- 
tive write and read pulses, At = 200 ns. 



FIG. 4: The experimentally derived ratio R T (ti,t2) [Eq. Q as 
a function of detection times (ti, £2) for the (1, 2) fields, with 
R T < 1 for classical fields. The left column gives R T (t\, £2) for 
At = 50 ns (top) and At = 200 ns (bottom), while the right 
column gives the associated statistical uncertainties. Bin size 
r = 30 ns. 



the write pulse, excitation is "stored" in the atomic en- 
semble until the readout. The production of correlated 
photon pairs should now be distributed along ti ~ 8t\2 
with width ~ Ati2- Instead, as shown in Fig. Hfc), 
p T (ti,<2) peaks towards the end of the write pulse (i.e., 
ti <; 100 ns), and near the beginning of the read pulse 
(i.e., 200 < t 2 < 300 ns). Early events for field 1 lead 
to fewer correlated events for field 2, as p T (ti,t2) decays 
rapidly beyond the line t^—tx = Td — 175 ns. The marked 
contrast between p T (ti, £2) for At = 50 and 200 ns results 
in a diminished ability for the conditional generation of 
single photons from excitation stored within the atomic 
ensemble [8J and, more generally, for the implementation 
of the DLCZ protocol for increasing At. The underlying 
mechanism is decoherence within the ensemble, as will 
be discussed. 

Fig. [2fb, d) displays q T (ti,t2) for independent trials 
j ^ k. q T (ti,t2) is expected to be proportional to the 
product of intensities of the fields 1 and 2, in reasonable 
correspondence to the form shown in Fig. EJb, d) for our 
roughly rectangular write, read pulses, but distinctively 
different from p T (ti, £2) in (a, c). 



To deduce R T {ti,t2) from Eq. ^ we next determine 
the joint detection probabilities p T (t\,ti) for field 1 and 
Pt (£2, £2) for field 2 from the same record of photoelectric 
events as for Fig. (a, c). Since the rate of coincidences 
for auto-correlations is roughly 10 2 times smaller than for 
cross-correlations for the (1,2) fields, we increase the bin 
size r to 30 ns to accumulate enough events to reduce the 
statistical errors to acceptable levels. Fig. shows the 
resulting time dependencies ofp T (ti,ti) andp T (i2,^2) for 
cases (/, 77). While the shape o£p T (ti, ti) associated with 
the write pulse does not change with At, the profile of 
Pr(t2,t2) from the read pulse is affected and exhibits a 
rise time that is ~ 3 times shorter for At = 200 ns than 
for At = 50 ns. This prompt rise in (b) is consistent 
with the observation that stored excitation is efficiently 
addressed at the beginning of the read pulse for non- 
overlapping write, read pulses, while the longer rise time 
in (a) results from continuous excitation and retrieval of 
atoms from the state \b) for the overlapping case. 

We employ the data in Figs. |21 El together with Eq. 
[flto construct the ratio R r {ti,t2), with the result pre- 
sented in Fig. 01 .13.]. Not unexpectedly, the trends 
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FIG. 5: Coherence time assessment. Experimentally acquired 
31,2 (black dots), theoretical description of the current exper- 
iment with K — 1.1 MHz (solid line), and the theoretical 
prediction for a spin polarized rriF = MOT (dotted line). 



for Rr{t\,t<2) closely resemble those of the joint proba- 
bility p T (ti,t2) for correlated pair generation previously 
discussed. As for the violation of the Cauchy-Schwarz 
inequality Rr(t\,t-2) < 1 for classical fields Il3i we 
observe maximal violations with i?™ ax = 292 ± 57 for 
At = 50 ns and Rf ax = 202 ± 60 for At = 200 ns 
(R T = 198 ± 33 in the neighboring bin). The relatively 
large errors in R T (t\,t2) arise predominantly from the 
uncertainties in p T {t\,ti) andp r (t 2 ,t 2 ) (Fig- El ,14]. 

The forms for p T (ti,t2) and Rr(ti,t2) for the cases 
At = 50 and 200 ns imply a decoherence process oper- 
ative on a time scale t& ~ 175 ns. To investigate this 
decay, we have performed a separate set of experiments 
with the delay At varied < At < 400 ns. For each 
At we determine the normalized correlation function g\ 2 
from the ratio of integrated coincidence counts to singles 
counts over the entire detection window (i.e., r = 200 
ns), with the results presented in Fig. [5] 

In Fig. [5J the initial growth of g\ 2 f° r small At is 
due to the finite time required to produce sequentially 
photons in the (1,2) fields, which is already evident in 
Fig. |3 More troublesome is the rapid decay of g\ 2 - 
A likely candidate responsible for this decay is Larmor 
precession among the various Zeeman states of the F = 
3, 4 hyperfine levels of the 65*1/2 ground level. 

To investigate t his p ossibility, we have extended the 
treatment of Ref. [l5j to include the process involving 
the read beam as well as the full set of Zeeman states for 
the F — 3, 4 hyperfine levels. The sample of Cs atoms is 
assumed to be initially unpolarized and distributed over 
the same range of magnetic fields as for the MOT. With 
write and read pulses that approximate those used in 



our experiment, we obtain an expression for the prob- 
ability p T {ti,t2) to generate a pair of photons at times 
(ti,t 2 ) for fields 1 and 2 as a function of the offset Ai. 
By summing contributions for all (ti,t 2 ) over the detec- 
tion windows, we arrive at the joint probability j>i, 2 (At) 
that we compare to the measured g\ t 2 (At) by way of a 
single overall scaling parameter for all At, as the rate 
of single counts in fields (1,2) is measured not to de- 
pend on At (to within 20%) The result is the solid curve 
in Fig. [SJ that evidently adequately describes the im- 
pact of Larmor precession on our experiment. The form 
of ^1,2 (At) strongly depends upon the inhomogeneity of 
Zeeman splitting across the MOT, which is described by 
the parameter K = ^b9f Lb/h, where L is the MOT 
diameter, b is the gradient of the magnetic field for the 
MOT, and gp is the Lande factor. The curve in Fig. [SJis 
the theoretical result for an initially unpolarized sample 
with K = 1.1 MHz as for our experiment (i.e., L « 3.6 
mm and b ss 8.4 G/cm). 

An obvious remedy for this dephasing is to eliminate 
the magnetic field altogether, as by transferring the sam- 
ple to a dipole- force trap [ll|. Alternatively, we have 
developed a scheme that should allow for long coherence 
times even in the presence of the quadrupole field of the 
MOT by utilizing only magnetic-field insensitive states. 
The write, read beams are polarized a± and are aligned 
along the z— axis of the MOT, which provides the quanti- 
zation axis. Atoms within the approximately cylindrical 
volume illuminated by these beams are initially spin po- 
larized into F = 3,m F = [l(|. The (1,2) fields are 
selected to be o± , which results in spin excitation stored 
in F — 4, m,F = 0. The prediction of our model for 
this new protocol for the same experimental conditions 
as before but now with an initially spin polarized sam- 
ple is shown as the dashed curve in Figure \5\ resulting 
in an increase of more than 3x in g\ 2 , and significantly 
extending the decoherence time to more than t^ ~ 10/is. 

In conclusion, we have reported the first observa- 
tions of the temporal dependence of the joint probability 
Pr{ti,t2) for the generation of correlated photon pairs 
from an atomic ensemble, which is critical for the pro- 
tocol of Ref. jj]. Our measurements of p r (ti,t 2 ) are 
an initial attempt to determine the structure of the un- 
derlying two-photon wavepacket 17]. The nonclassical 
character of the emitted (1,2) fields has been tracked 
by way of time dependence of the ratio -R T (ii,t 2 ), with 
i? max = 292 ±57 ^ 1. Decoherence due to Larmor preces- 
sion has been characterized and identified as a principal 
limitation of the current experiment. A new scheme for 
effectively eliminating this decay process has been pro- 
posed and analyzed, and could be important for the ex- 
perimental realization of scalable quantum networks [4| 
as well as for an improved source for single photons y|. 
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